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Distortions of the oxygen sublattice couple to e g orbitals of 
Mn 3+ and drive a cooperative Jahn- Teller (orbital ordering) 
transition in LaMnOs. A simple model for this transition 
is studied. Without further adjustment, the model predicts 
the shape and stability of small (anti- Jahn- Teller) polarons 
which form when holes are doped into the material. This 
leads to a new description of the lightly doped insulator, the 
anti-ferromagnetic to ferromagnetic transition, and the metal- 
insulator transition. 

71.38.+i,71.30.+h,75.30.Vn 



I. INTRODUCTION 

The doped manganites Ri-^S^MnOs (where R is typ- 
ically La and S is typically Sr or Ca) have a fascinating 
(T,x) (temperature, concentration) phase diagram, in- 
cluding "colossal magnetoresistance" JjJ when T 250K 
and x «0.20. Experiment (transport optical ||, 
diffraction §|,EXAFS §], isotope studies @) and the- 
ory indicate the appearance of polarons at this point 
in {T,x). It is less well appreciated that small polarons 
are essential to explain the insulating phase at smaller 
x and 0K< T <750K §. We study a simple model for 
the cooperative Jahn- Teller (JT) transition, and use the 
model to predict properties of small polarons in lightly- 
doped material, including how they affect magnetic order 
and the metal-insulator transition. By working in the 
limit of large on-site Coulomb repulsion U, we find that 
properties of the polaron are simple enough to describe 
analytically, with small perturbative corrections. 

There is disagreement over the relative role of 
Coulomb, magnetic, and electron-phonon effects. We 
offer a simple unified picture in which the relevant en- 
ergy scales in descending order are: (1) Coulomb inter- 
actions are inactive after establishing the dominant Hub- 
bard and Hund energy scales; (2) electron-phonon inter- 
actions drive orbital ordering by the JT mechanism; (3) 
when doped, electron-phonon interactions via orbital or- 
dering give small anti-JT polarons; and finally (4) orbital 
ordering disrupted by polarons is responsible for the un- 
usual magnetic phases (described in Sec. III.) 

A canonical "Jahn- Teller polaron" Q is the excess elec- 
tron in BaTi03 [[lO) , which sits in a triply degenerate tig 
level. A local distortion of the lattice splits the degener- 
ate levels. Beyond a critical coupling strength, this low- 
ers the energy and traps the electron in a small polaron 
state. 
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FIG. 1. Base plane (x-y plane) of Jahn- Teller distorted 
LaMnOs. Rotations of oxygen octahedra are omitted, and 
distortions are exaggerated. Oxygens are displaced only along 
Mn-O-Mn bonds. 

By contrast, the polarons in LaMnOs are "anti-JT po- 
larons." When pure, the valence is Mn 3+ , with d A con- 
figuration in the high-spin state t^g^e 1 ^ with the doubly- 
degenerate e g level singly occupied. Below the JT struc- 
tural transition at Tjt=750K [jll), oxygen octahedra dis- 
tort as illustrated schematically in Fig. [jj lifting the 
e g degeneracy and lowering the energy of the occupied 
orbitals. When a hole is added, we show that a small 
polaron is formed by locally "undoing" B the JT distor- 
tion, pinning the hole onto a Mn 4+ site with a filled tig 
shell. 



II. MODEL HAMILTONIAN 



The Hamiltonian we use is H = TL t + H cp + Hl + Ti.v- 
The first term represents hopping of Mn e g electrons to 
nearest neighbors. A simple way to derive this term is to 
introduce an overcomplete basis 



ip x = 3x - r , vb y = 3y - r , ip z = 3z 



(1) 



each pointing toward the two nearest Mn neighbors along 
one of the Cartesian axes. Note that ip x + tp y + ip z = 0. 
In the two-dimensional e g space, these basis vectors lie 
at 120° to each other, as shown in Fig. |[ The usual 
orthogonal basis, ip2 — {ipx — Vv) and tp3 — ip z , is 
more complicated because there is no element of the cubic 
rotation group which transforms ip2 into ^3 . Then for the 
hopping Hamiltonian we choose: 



Ht = t 



L± 



{[ct(£)c x (i±x)] + [x^y} + [y^z]} (2) 



1 



where t numbers manganese sites, and I + x numbers 
the Mn neighbor to the right. After re-expressing Ti t in 
terms of the orthogonal orbitals ip2 and ^3, we recover the 
correct nearest-neighbor two-center Slater-Koster 
model, with overlap integrals t = (dda), (ddir) not en- 
tering due to symmetry, and (ddd) = 0. 




FIG. 2. The two-dimensional space of e g orbitals, with ver- 
tical and horizontal axes being the usual orthogonal basis 
functions 1^2 and ipa. The symmetrical overcomplete basis 
functions ip x , ip y , and ip z , Eqn. are shown as dashed lines 
and the symmetrical orthogonal basis ipx and ipy, Eqn. 
are shown lying at 45° in the second and third quadrants 



The only lattice degree of freedom is oxygen motion 
along the direction of the bonds to the nearest two Mn 
atoms, with a harmonic restoring force —Ku x . Variables 
Qx{£) = u x {£ + if) — u x (£ — if) measure the local ex- 
pansion of oxygens on the x-axis around the £-th Mn 
atom. This expansion lowers the energy e x of orbital ip x 
by de x /dQ x — — 4g/\^3. We work in adiabatic approxi- 
mation (oxygen mass =00). This gives 



He 



^§ E {[4(*M*)Q*(*)] + [x - y] + [y - z]} 



(3) 



same model, but with longer range forces, was used by 
Millis 0] for LaMn0 3 . 

To analyze the solution, it is convenient to make a 
45° rotation in (tp2, ^3)-space, to a new orthogonal basis, 
shown by the solid arrows in Fig. || 

<l>x= -^=(V2-^3) = ^=((V3 + l)Vx+(V3-l)Vy) 

</V = -4=fofc» +V*) - ^= ((V3- l) ^ + (V3 + l) ip y ) (6) 

The orbitals ipx and ipy point strongly in the x and y 
directions, and are equivalent under a 90° rotation in real 
space, as shown in Fig. [I]. 

In the new basis the electron-phonon term in the 
Hamiltonian (0) is conveniently split into Jahn- Teller and 
breathing parts, H cp — Hjt + Tibx- 



Hjt 



x 



Q2W 

QaV) 



-<M*) 



exit) 
c Y (l) 

(7) 



I H hl = -V2(l + P)gY,Qi(^x{^x(t) +Ml)c Y (l)] 



(8) 

where the breathing amplitude is Q\ = y/2/3(Q x + Q y + 

Qz), Qi is Q x ~ Q y , and Q 3 is (2Q Z - Q x - Q y )/V3, 
in standard Van Vleck [jl5[ notation. A non-zero f3 was 
introduced by Millis, to represent additional charge cou- 
pling to the breathing mode. To simplify the model, we 
use = 0, which makes Eqns. identical to Eqn. [| 

The hopping Hamiltonian in the new basis is: 



H t = t 



T, 




£,5=x,y,z 



,T,= 




1 1 

! I 



(4) 

When two orbitals are occupied on one site, one has to 
pay Coulomb energy U : 



(5) 



We take into account t2 g spins, by defining the c 1 opera- 
tors to create electrons whose spins (because of the large 
Hund's rule energy) are parallel to the S=3/2 core spin. 
Hopping (Eqn. |J) then operates only between adjacent 
Mn sites with parallel spins. 

This model generalizes to e g electrons the model of 
Rice and Sneddon |Tj| for s electrons in BaBiC^. The 



III. GROUND STATE SOLUTION 

We have solved this Hamiltonian for zero doping [x = 
0) in two opposite limits: U/t small, by a Hartree-Fock 
calculation, and U jt infinite. The two limits give simi- 
lar answers. The latter case seems to us more realistic 
and has another advantage: since hopping is prevented 
by unit occupancy of all sites, the solution does not de- 
pend on the magnetic order. We postpone discussing the 
Hartree-Fock solution until the end of this section. 

A distortion (Q2,Q3) = Q(cos0, sin#) exp(i<f • £) is in- 
troduced. We minimize elastic and electron energy to 
find optimal distortion parameters Q, 9, and q. The op- 
timal distortion has wavevector q — (tt, n, n). Unless we 
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add anharmonic or strain terms to select a direction in 
(Q2, Q3)-space, the energy is independent of 9. Exper- 
iment shows that the actual ordering is Q2-type, so we 
simply choose this distortion, and avoid having an extra 
term in the Hamiltonian. 

The ground electronic state which corresponds to per- 
fect Q2-type orbital order is 

A B 

iJT>=n c xwn4-( f )i o> < (10) 

i £' 

where A and B label sublattices where the phase of the 
orbital order exp(iq ■ f) is ±1 respectively. The energy 
< Y£\H\ JT >= -NgQ + NKQ 2 /16 has minimum value 
E/N = -4Tt at g = Bg/K as shown on Fig. |. This 
JT phase is insulating; the gap to charge excitations is 
approximately U ~ 6eV. Electron-phonon effects are con- 
veniently expressed in terms of the dimensionless param- 
eter r = g 2 /Kt which we estimate to be «0. 25-0. 35. 
Orbital excitations |l6| require only energy 16Fi «2 eV, 
and spin excitations occur down to low energies, with 
energies «50K determined by balancing terms of order 

As explained by Goodenough jll| , the Q2-type orbital 
order (observed below Tjt=750K) gives a layered struc- 
ture (shown in Fig. |l|) which in turn fixes the spin order 
which sets in below the Neel temperature T/v =140K. 
The magnetic structure seen |fl7| at small x is "antiferro- 
magnetic A" (AFA), with spins aligned ferromagnetically 
within the layers and antifcrromagnctically perpendicu- 
lar. The source of the ferromagnetic in-plane exchange 
(JiS 2 =3.32 meV E^]) is orbital order which favors vir- 
tual hops from filled A sublattice ipx orbitals to empty B 
sublattice tpx with spin parallel to avoid a Hund penalty. 
When orbitals are not ordered, the Hund penalty is out- 
weighed by the greater multiplicity of hops which tig- 
electrons can make onto empty ^-states on antiparallcl 
neighbors. The antiferromagnetic c-axis exchange with 
J2S' 2 =-2.32 meV is explained this way. 

It turns out that none of the energetics described so far 
(for Q2-type orbital order) depend on the z-component of 
the orbital ordering wave-vector q, which experimentally 
is (ir, 7t, 0) rather than (it, it, it) as would be preferred for 
any orbital order other than pure Qi. Again, we simply 
adopt this order, without specifying the additional term 
in the Hamiltonian. 

In the opposite limit U ft small, hopping energy is dom- 
inant, and the results depend on the magnetic order. We 
consider the two cases of ferromagnetic (Ferro) and AFA 
order. In latter case we turn off hopping in the z direc- 
tion (T z = in Eqn. fjl) A gain, a distortion (Q2,Qs) = 
(2(cos#, sin#) exp(i<f • £) is introduced, and we minimize 
elastic and electron energy in Hartree-Fock (HF) approx- 
imation. At half-filling (x = 0) of the e s -bands, there 
is perfect nesting at wavevector q — (vr, tt, if), because 
of a symmetry of the bands of the undistorted crystal, 
ei(fc) = — £2(fc + q). This fixes the optimal distortion 



at q = (%, 7r, 7r). The resulting HF energy is almost in- 
dependent of the angle 6 in (Q2,Q3)-space, favoring the 
(^-distortion by one part in 10 4 of the JT energy. The 
calculations are shown in Fig. 0, which also shows the 
critical value of U/t at which the weak coupling (HF) so- 
lution and the strong-coupling (U /t — > 00) solutions have 
equal total energy. For LaMn03, we estimate U/t w 12. 
The lower panel of Fig. || shows that for this choice, the 
U — > 00 approximation is better. The upper panel shows 
that the JT distortion does not depend much on U for 
r > 0.25. 
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FIG. 3. Critical value of on-site repulsion U which sepa- 
rates weakly from strongly-correlated solutions. The upper 
panel shows the Jahn- Teller distortion p — gu/t as a function 
of r = g 2 /Kt for weak and strong-coupling solutions. 



IV. SMALL POLARON 

What happens when a hole is added? It can go into 
any of the states \£A > and \£'B > which are occupied 
in Eqn. When T = the hole is free to hop among 
these states, no matter how large is U, since spins are 
aligned in the planes. When T > it costs energy gQ = 
8Tt in lost JT energy to remove the electron. There are 
two different ways to regain some energy, (a) The hole 
can delocalize, forming a conducting Bloch state, with no 
relaxation of the oxygen coordinates Q(£). For U = 00 
the wavefunction in first approximation is 

A B 

^k=c A ^ e iU \£A > +c B ^ \?B > (11) 

e £' 

where the N states \£A >, \£'B > are obtained from 
I JT > by putting the hole on a single site; \£A >= 
c x {£)\JT > and \£'B >= cy(£')\JT >. The resulting 
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hopping energy is ek/t — (l/2)(cos k x + cosk y ) plus an 
additional term — cos k z in the Ferro case due to hopping 
in z direction (prohibited in AFA case.) The minimum 
energy of the extended hole state is -Eh.cxtA = 8r — 1 
(AFA) and 8r - 2 (Ferro). (b) If the distortions Q e 
are locally readjusted, a bound state can be formed, 
which will be non-metallic because a small impurity 
potential will pin it. In first approximation, put the 
hole at a single A site (neglecting hopping for now.) 
The nearest i-oxygens, instead of being displaced out- 
wards by u — 2g / K , should now displace inwards by 
u = — (a/4/3— l)g/K; the y-oxygens, formerly displaced 
inwards by u = —2g/K, now displace slightly further 
inwards, by u — — (a/4/3 + l)g/K; the z-neighbors, for- 
merly undisplaced, now displace in by u — —^A/Sg/K. 
The pattern of this anti-JT polaron is shown in Fig. ||. 
The energy i?h, P oiA needed to make the hole is reduced 
from 8r in lost JT energy to 2r. This energy comes 
partly from a reduction in the strain cost and partly from 
energy of breathing. Setting i^h.oxt = -Eh,poi, we con- 
clude that small polarons are stable for T > 1/6 in AFA 
phase, and that ferromagnetism, by enhancing the der- 
ealization energy, prevents small polarons until T > 1/3. 
These energies are shown as the solid curves in Fig. EL 
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Electron-Phonon Coupling r 
FIG. 4. Lowest hole energy in the AFA state of pure 
LaMnC>3 as a function of coupling constant V, for the trial 
delocalized and localized solutions. Solid curves are lowest 
order variational solutions and dashed curves include pertur- 
bative corrections. The dash-dot curve is an improved pertur- 
bative solution (exact subspace diagonalization) for the local- 
ized case. The crossover points at each level of approximation 
are shown in boxes. The inset shows (in exaggerated form) 
the shape of anti-Jahn- Teller polaron (dashed lines) in the 
x-y plane. The relative positions of O atoms (at vertices of 
rhombus) in the pure JT (solid line) state and in the polaron 
state are accurate. 

The three largest errors in the calculation of the hole 
energy are (1) the localized hole can spread somewhat 
onto neighboring Mn atoms; (2) the Hilbert space for 
both localized and delocalized hole wavefunctions should 
include states with orbital defects (sites occupied singly, 
but by the orbital of wrong orientation); (3) a hole in 
AFA phase will lower its energy further by causing spin 



canting on Mn atoms in adjacent planes, thus permit- 
ting intcrplanar hopping, and eventually driving an AFA 
to Ferro transition. We have calculated these effects by 
perturbation theory. Details for processes (1) and (2) 
are given in the Appendix, and process (3) is further dis- 
cussed in the next section. 

The probability that a small hole polaron will be found 
away from the central site is a/(1.0366 + a) where a = 
0.0059/r 2 (AFA) and 0.0099/r 2 (Ferro). At the smallest 
r where small polarons are stable, the probability is 19% 
and 10%. This shows that the small polaron is very well 
localized. The first two effects cause the energy to be 
lowered by CtT + C't/T with C=0 for delocalized states 
and 0.49 for localized states, and C"=(0.11,0.14) in (AFA, 
Ferro) for delocalized states and (0.11, 0.17) for localized. 
The shifts are small for values of T > 0.2. The critical 
value of r for small polaron formation changes from r c 
I fi to 0.157 (AFA) and from 1/3 to 0.294 (Ferro). The 
corrected energies are shown as the dashed curves in Fig. 

I- 

V. MAGNETIC TRANSITION 

Effect (3), the spin-canting problem, has usually been 
discussed on the assumption that holes are delocalized 
p9[ . Treating spins classically, the additional dereal- 
ization energy gained when adjacent layers realign from 
180° to 180° - 9 is t sin(0/2) per hole; the exchange cost 
of realignment is J2S 2 (1 — cos(8)) per Mn atom. The 
optimum canting angle is sin(6>/2) = x</4| ^l^ 2 , which 
gives a small critical concentration x c — A\J2\S 2 /t for 
complete rotation to the Ferro phase. Localized holes 
tilt spins on only near-neighbor Mn atoms in adjacent 
planes, gaining less hopping energy than delocalized holes 
because the electron hops into an anti-JT-distorted hole 
site. The energy gain is (39/640)(t/r) sin(0/2). If we 
neglect the rotation of any except first-neighbor spins, 
the magnetic energy loss around the localized hole is 
(8Ji + 2| J 2 \)S 2 (1 — cos(0)) per hole, giving the optimum 
local rotation sin(0/2) = 0.49(eV" 1 )t/r. For t > 2MT 
eV, the adjacent spins are completely flipped. Compar- 
ing the magnetic energy loss per spin with the energy loss 
of the Ferro state, we find a critical concentration plotted 
in Fig. ^. Experimentally the AFA/Ferro phase bound- 
ary occurs at x c =0.08 for Sr-doped LaMnC-3 Q and at 
x c =0.15 for Ca doping |l7]]. These values are easily rec- 
onciled with the localized hole picture, but demand too 
small a value of t in the delocalized picture. The smaller 
Ca ion causes larger Mn-O-Mn bond angles which results 
in a smaller t and a larger x c . 
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FIG. 5. Critical concentration x c for the AFA to Ferro 
transition as a function of coupling constants t/Y. 

VI. METAL-INSULATOR TRANSITION 

At doping concentrations x approaching x c =0.20, 
LaMnOs becomes metallic. A possible analog system 
with simpler but similar physics is doped BaBi03. When 
Ba atoms are replaced by K atoms, the system remains 
insulating up to a concentration :r=0.40. The most likely 
and successful explanation for the insulating behavior is 
formation of bi-polarons. The pure material has Bi atoms 
in the nominal Bi 4+ valence state. Chemically Bi prefers 
valences Bi 3+ and Bi 5+ . The pure material has alternat- 
ing oxygen displacements inward or outward in the Qi 
breathing pattern. This stabilizes charge ordering which 
is interpretable as alternation of Bi 3+ and Bi 5+ . When 
K atoms replace Ba, extra electrons are released from Bi 
atoms, creating holes, and converting some of the Bi 3+ 
ions into Bi 5+ . The new Bi 5+ ions created by doping are 
small hole bi-polarons, and are stabilized by relaxation of 
the Qi coordinates of oxygens. A microscopic description 
has been given [^jj using the Rice-Sneddon Hamiltonian 
Jk|, which is just a simplified version of the Hamiltonian 
used here, with a single Bi s-band in place of the two 
Mn e g -bands. The Hubbard U is much smaller for Bi s 
orbitals than for Mn e a 's, so bi-polarons are the stable 
defect of the Peierls-type charge order parameter, while 
single polarons are the stable defect of the JT orbital 
order parameter of LaMnC>3. 

The metal-insulator transition in BaBiC>3 was ad- 
dressed in a remarkable calculation by Yu, Chen, and 
Su (verified by Kostur and Allen p3[). In their cal- 
culations, bipolaron defects form spontaneously in the 
Peierls order parameter up to very large concentrations, 
but when the concentration is too high, the defects de- 
stroy the Peierls order completely, and the ground state 
switches from defective dimerized insulator to undis- 
torted metal. In the work of Kostur and Allen |23| it is 
observed that the defective dimerized state (a "bipolaron 
glass" phase) is probably an Anderson insulator. The 
bipolarons are dense enough that localized single-particle 
states fill up the Peierls gap, and an Efros-Shklovskii 
pseudogap p3|] rather than a clean gap separates the filled 



and empty states at the Fermi level. The insulating state 
might have a non-zero linear specific heat coefficient 7. 
However, the metal-insulator transition is not an Ander- 
son transition - the localized states do not persist away 
from the Fermi level in the metallic state, but disappear 
altogether. 

The transition in LaMnOs seems to us to be similar, 
with polarons instead of bipolarons, and JT rather than 
Peierls order being destroyed at the phase transition. We 
have not attempted a calculation, which would be more 
difficult for LaMnOs than for BaBi03 because of the 
large value of U. Schematically the energy of JT glass 
phase and metallic phase would be something like 

E 3T /Nt = -[4- a{U)x - b{U)x 2 + ■ ■ -]T 
E meta .i/Nt = - [c(U)x + d(U)x 2 + ■■■] (12) 

where coefficients a, 6, c, d, etc. are unknown functions 
which depend strongly on U and also on the magnetic 
state. Our calculations above give the value of a at U = 
00 to be 2 plus small corrections. If we could neglect 
bx relative to a and dx relative to c, then the critical 
concentration would be x c = 4r/(ar + c). The observed 
x c « 0.2 then requires that the hopping energy gain — ct 
in the correlated metal phase must be of order -5t, using 
r = 0.3, whereas we find that at U = the value of c in 
the ferromagnetic phase is only 3. This shows that the 
polaron repulsion b plays a significant role. 

VII. DISCUSSION 

Numerous effects are left out of the model. Macro- 
scopic strain, tilting of oxygen octahedra, and additional 
zero-point magnetic and non-adiabatic lattice fluctua- 
tions all could be added and would affect the T = 
properties discussed here, but we believe that the ef- 
fects we did include are the most important ones. A 
T > treatment of all these effects is a bigger challenge. 
It would also be interesting to extend the T = cal- 
culations to higher doping where a remarkable variety 
of textured phases is being unraveled p5[ . The proper- 
ties of LaMnOs-related materials are incredibly rich, yet 
surprisingly understandable, in contrast to certain other 
transition-metal oxide systems. 
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APPENDIX: PERTURB ATIVE CORRECTIONS 

In Sec. IV the energy of the hole in JT-distorted 
LaMnC>3 is computed in lowest order, for the two cases 
of localized and delocalized hole states. Here the method 
used for perturbative corrections is explained. First we 
discuss the localized case. 

Denoting the pure crystal ground state Eqn. as 
| JT >, the strictly localized hole on an A site at the origin 
is cx(0)|JT >. For clarity we use both a shorthand and 
a pictorial labeling scheme: 



\h0 >=c x (0)|JT >= 



I 

I I 
I 



(Al) 



where hO means "hole at the origin" . This is the zeroth 
order (unperturbed) state of the localized hole. As shown 
in section IV, relaxation of oxygen positions reduces the 
energy e(hO) =< hO\H to t\hO > from 8Tt to 2Tt. Our aim 
is to calculate the two corrections, 



5Ei 



E 

i=l 



\(i\HjT\hO)\ 2 
eQiO) - e(i) 



(A2) 



(A3) 



Included in the sums are all basis functions which cou- 
ple to \h0 > by either H jt or Tit- All our basis states 
have the same relaxed oxygen positions optimized for the 
\h0 > state. 

The JT operator Tt jt couples \h0 > to six states \i > 
with orbital defects on nearest neighbor Mn atoms. The 
states \i > come in three types. First consider 



hO,ox >= c^(i)cr(i)c x (0)|JT >= 



I - 



(A4) 



and the related state \h0,o — x >. The label \hs, os' > 
has the meaning "hole at site s, orbital defect at site 
I'." When there is no orbital defect, the second part of 
the label is omitted. These states have energy e(hO, o ± 
x)-e(hO) = (14 - 4/V3) Tt, reduced from the value 16rt 
of a distant orbital defect by relaxation of the displaced 
oxygen lying between the origin and the ±i Mn. These 
states are coupled to | hO > in first order by the TL jt term 
proportional to Q 3 (±x) — — (2/3 + l/Vo) g/K. Second 
are the states \h0, o±y > with the hole at the origin and 
the orbital defect at the ±y sites. These have a slightly 
increased energy denominator (14 + 4/Vo) Tt, and cou- 
pling amplitude to the \h0 > state caused by the oxygen 
relaxation Q 3 (±y) = - (2/3- 1/V3) g/K. Third, the 
states \h0, o ± z > have the same energy denominator 
16Tt as a distant orbiton, because the only neighboring 



oxygen with an altered position lies along the z axis, 
generating the term Q 3 (±z) = 4:g/3K which causes off- 
diagonal coupling. Adding the three types of correction, 
we get the value for the first perturbative correction to 
the localized hole energy, 



628 
1287" 



Tt 



-0.488rt. 



(A5) 



This correction is independent of magnetic ordering. 

The hopping Hamiltonian TL t couples state \h0 > to 
twelve states, which come in two categories of six each. 
One example of each category is: 




\hx >= c Y (x)\JT >- 



\hx,oO >= c f Y (0)c x (0)c Y {x)\JT >= 



I 

I I 
I 



(A6) 



(A7) 



The first six have the hole displaced to a first neighbor 
site, with the origin occupied by an electron in the "cor- 
rect" orbital; the second six have a misoriented orbital at 
the origin. The hops without orbital defect in the ±(s, y) 
directions cost energy e(j) — e(hO) = 37ri/3, while the 
hops in the ±z direction (forbidden in the AFA magnetic 
state) cost 40rt/3. Moving a hole to a remote site would 
cost slightly less energy 12Tt. Hops which leave an or- 
bital defect at the origin cost an extra 8Tt. The coupling 
magnitude is determined by the off-diagonal part of the 
matrix of Eqn. ^| for hopping without creating orbital 
defect at the origin and by the diagonal parts otherwise. 
The net result is 



^(Ferro) 
SE 2 (AFA) 



3 


3 


21 


148 H 


-80 + 


244 


3 


21 > * 




148 n 


" 244 j 


If- 



3 

128 
-0.106- 



-0.167^ 

(A8) 



These calculated shifts are shown in the AFA case as 
a dashed line in Fig. ^. The shifts are similar in magni- 
tude to the unperturbed energy. Therefore we repeated 
the calculation with an exact diagonalization in the same 
subset of 18 states which couple to \h0 >. The answer, 
shown as the dot-dashed line in Fig. ||[ is not greatly 
different from the perturbative answer in the relevant 
regime of parameters. 

In section IV the case of the delocalized hole was 
treated by staying within the subspace of states c(£) | JT > 
where c(£) removes an occupied JT orbital, and then di- 
agonalizing the off-diagonal effects of 7i t by Fourier trans- 
formation to Bloch states. The hole is then put into the 
lowest energy Bloch state. This calculation has omitted 
effects caused by the fact that 7i t also allows hops which 
leave behind a single orbital defect. We now correct for 
this perturbatively. Our Hilbert space has two subspaces. 
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The Bloch states lie in the iV-dimensional subspace with 
a single hole and no orbital defect. We must add a 6N- 
dimensional space in which the hole has an orbital defect 
on an adjacent atom. These two subspaces are coupled 
by Ti t . The Schrodinger equation has the structure 



Shrekala, R. Rames h, M. Rajeswari, and T. Venkatesan, 



{ K 



e n' t 



011 



0. 



(A9) 



Each element of the ipi subspace (no orbital defect) is 
coupled to six elements of the ipu subspace by hopping 
terms. The prime on H.' t is used to designate the part 
of Ht which creates an orbital defect. The JT energy is 
constant (8Tt hole creation energy) in the ipi subspace 
and higher by the JT gap (16Tt) in the ipu subspace. In 
our perturbative treatment we leave out the off-diagonal 
effects of H.t interior to the ipu subspace. The problem 
is then equivalent to an effective Hamiltonian 



Wcff = Kjt + T~Lt — Wj(Wjt - E) 1 TL' t 



(A10) 



in the subspace. Second-order perturbation theory 
uses the JT gap IQTt as the en ergy denominator (Hjt — 
E). Then the last term in Eqn. AlC just gives a constant 
shift on the diagonal. The value depends on the magnetic 
state. 



-Q.141- 



6E 3 (AFA) = - 



64T 



-0.109- 



t 



This shift is shown as a dashed line (in the AFA case) in 
Fig. § 
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